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Abstract 



We study a model for an argon-like fluid parameterised in terms of a 
hard-core repulsion and a two- Yukawa potential. The liquid-gas phase be- 
haviour of the model is obtained from the thermodynamically self-consistent 
Ornstein-Zernike approximation (SCOZA) of H0ye and Stell, the solution of 
which lends itself particularly well to a pair potential of this form. The pre- 
dictions for the critical point and the coexistence curve are compared to new 
high resolution simulation data and to other liquid-state theories, including 
the hierarchical reference theory (HRT) of Parola and Reatto. Both SCOZA 
and HRT deliver results that are considerably more accurate than standard 
integral-equation approaches. Among the versions of SCOZA considered, the 
one yielding the best agreement with simulation successfully predicts the crit- 
ical point parameters to within 1%. 
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I. INTRODUCTION 



Much attention has been paid in recent years to the hard core Yukawa (HCY) potential 
as a model for the pair interactions of fluids [0J . Interest in the potential is motivated on the 
one hand by its relevance to solvent averaged interactions in polyelectrolytic and colloidal 
particles and, on the other hand, by its analytical tractability in the context of liquid state 
theories such as the Mean Spherical Approximation (MSA) and the Self Consistent Ornstein- 
Zernike Approximation (SCOZA). Recent studies of the HCY fluid can be found in 0|| and 
references therein. 

For simple fluids interacting via dispersion forces, however, the bare HCY potential fails 
to provide a realistic representation of the interactions, which are much better modelled 
by the Lennard- Jones (LJ) potential. Unfortunately owing to its mathematical structure, 
the latter is less amenable to direct study by the MSA and SCOZA than is the HCY 
potential. Nevertheless, it has long been appreciated that a LJ-like simple-fluid potential 
can be replaced by a potential with a hard core plus a linear combination of two Yukawa tails 
with no appreciable loss of agreement between the experimental and model equations 
of state. Such a potential permits the theoretical study of simple fluids whilst retaining the 
convenient mathematical properties of the HCY potential. 

Several different representations of a LJ-like fluid in terms of a hard-core plus two Yukawa 
fluid (HC2YF) have been proposed in the literature . The one we use here is very similar 
to that recently given by Kalyuzhnyi and Cummings |J: 

oo r < a , 

exp[— z\[r — o)\ exp|— z 2 (r — a)\ r > a . 

Here a and e set the length and energy scales of the model and coincide with the zero and 
the well depth of the LJ potential respectively. The parameters At = 1.6438a, z\ = 14.7cr _1 , 
A 2 = 2.03a and z 2 = 2.69er~ 1 are chosen to fit the LJ potential and match the second virial 
coefficient to that of the LJ fluid. Eq. ( |1 . 1| ) differs slightly from the original parameterisa- 
tion by Kalyuzhnyi and Cummings insomuch as these authors determined the hard-sphere 
diameter <ths as a function of temperature via the Barker-Henderson procedure [|1(J in or- 
der to represent the soft-core repulsive interaction which results from the Weeks- Chandler- 
Henderson separation of the LJ potential [T1J. We instead hold the hard-core diameter fixed 
by setting o"hs = cr for conceptual simplicity and clarity in making comparison with our 
simulation results. A comparison of Eq. ( |1 . 1| ) with the LJ potential is shown in Fig. [I]. Ka- 
lyuzhnyi and Cummings have used the MSA, Percus-Yevick approximation and reference 
hypernetted chain theories to obtain the liquid-gas phase behaviour of the HC2YF specified 
by their parameterisation. Their results showed good agreement with existing literature 
data for the LJ fluid. 

In the present work, we extend the work of Kalyuzhnyi and Cummings || by determining 
the liquid-gas phase behaviour of the HC2YF using the SCOZA of H0ye and Stell. Unlike 
those considered in ref. ||, this theory is neither mean-field nor mean-spherical like and on 
the basis of a similar study already performed for the HCY fluid, it can be expected to give a 
superior account of fluid properties in the near-critical regime. We compare the predictions 
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of the SCOZA for the the critical point and the coexistence curve to those obtained by other 
approaches, including the renormalization-group based hierarchical reference theory (HRT) 



of Parola and Reatto [T2], and to new high resolution simulation results for the HC2YF. 
We find that one version of the SCOZA considered provides a remarkably accurate critical 
point and coexistence curve. The critical density and temperature predicted by the theory 
agree with the simulation results to within 1%. 



II. THEORY 

The SCOZA deals with two-body potentials which, like that of Eq. ( |1.1|) , consist of 
a singular hard-sphere repulsion with diameter a and a longer-ranged tail w(r). As is 
customary in integral-equation theories, this approach is based upon the Ornstein-Zernike 
(OZ) equation linking the two-body radial distribution function g{r) to the direct correlation 
function c(r). A closed theory is obtained by supplementing the OZ equation with an 
approximate relation involving g(r) and c(r). In its simplest form, the SCOZA amounts to 
setting: 

g{r) =0 r < 1 , 

(2.1) 

c(r) = K(p, (3)w{r) r > 1 , 

where p is the number density of the system, (3 = 1/ (k^T) is the inverse temperature, and the 
hard-sphere diameter has been set equal to one. This closure resembles that adopted in the 
MSA, except that the amplitude K of the direct correlation function outside the repulsive 
core is regarded as an unknown state-dependent quantity, to be determined in such a way 
that consistency between the compressibility and the energy route to thermodynamics is 
enforced. This constraint amounts to requiring that the reduced compressibility Xred and 
the excess internal energy per unit volume u satisfy the condition: 

=<V- < 2 - 2 > 

where it is understood that Xred is obtained by the compressibility sum rule as the structure 
factor S(k) evaluated at k = 0, while u is obtained by the energy equation as the spatial 
integral of the tail interaction w(r) weighted by the radial distribution function g(r). If the 
closure (|2.1| ) is used for the correlations, the consistency condition ( |2.2j ) yields a closed partial 
differential equation (PDE) for the function K(p,j3). In the present case, implementing this 
scheme is made simpler by taking advantage of the analytical results obtained for the OZ 
equation with the closure fl2.1|) when the tail potential w(r) has a two- Yukawa form like 
that of Eq. ( |1.1| ). These enable one to obtain Xrcd as a function of p and u. The function 
Xred(p, u) can then be used in Eq. ( |2.2|) by taking u instead of K as the unknown quantity. 
The algebraic manipulations are similar (although not identical) to those we performed in 
Ref. 0, and are based on the results obtained in Ref. |T3| . We start by introducing the 
quantity: 

/ = (l-0^, (2-3) 

V Xred 
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where £ = np/6 is the packing fraction. This can be written as: 



./ 



[z\ ~z 2 ) + ( 72 - 71) - z| 

4 [(zi/z 2 ) 72 - (^Ai) 7i] z x z 2 [{z x /z 2 ) 7 2 - (^/zO 71] 



7172(72 - 7i) 



2 • 



(2.4) 



where q — (1 + 2£) 2 /(l — £) 2 , and 71, 72 are state-dependent functions. By using Eqs. 
and (|2lD , Eq. Q becomes: 



2/ 



s 



'9/ 
972 



P>7i 



du 



' d 2 u 



(2.5) 



where the partial derivatives of / with respect to 71, 72 are straightforwardly obtained from 
Eq. (J2.4Q. In order to obtain a closed PDE for u, one then needs to express 71, 72 as a 



function of p and u. To this end, we recall the definitions of 71, 72 |13 

2- Vq 



A + 2zj- z 2 Tjlj-l 
2(2 + Zi ) aJi-1 



1,2), 



(2.6) 



where the quantities are known and depend only on the inverse range Z{, and Ii is 

given by: 



/+00 
dr rexp[— Zi (r — 1)] g(r) (i — 1,2). 



(2.7) 



For the two- Yukawa tail potential of Eq. ( |1 . 1| ) , the integrals Ji, J2 are straightforwardly 
related to the internal energy per unit volume u: 



u 



pe {A x h - A 2 I 2 ) ■ 



(2.8) 



If Eq. ( |2.6D is used to express Ii as a function of 7$ and the result is substituted into Eq. 
one obtains an expression for u as a function of p, 71, 72. This can be inverted algebrically 
to obtain 72 as a function of p, 71, and u. One is then left with the task of expressing 71 in 
terms of p and -u. If we set: 



47i 

z 2 

y = Vq--r- 

472 



we find by applying the results of Ref. [13| that the following equation hold: 



A 2 zla 2 l (2 + z 1 ) 2 *{2-y/q- ai )(y/q 



x — z\ 



{4(z 2 2 -4y 2 )(y-x) 



z\ + 4(y 



2 2 



)]} 



- A x z\ o\{2 + z 2 ) 2 [4(2 -^q-a 2 ){^/q-y)- z 2 2 ] {a{z{ - Ax 2 ){y - xf 

4(y 2 -x 2 )]} = 0, 



Z 1 2^9 



(2.9) 
(2.10) 



(2.11) 
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where i = 1,2, are known functions of the inverse-range parameters Zj. Note that this 
equation does not contain the unknown amplitude K(p, (3). By writing x and y as functions 
°f 7i> 72 via Eqs. (|2.9|), ( |2.10| ), and subsequentely 72 as a function of p, 71, w, we finally 
obtain an equation for 7! as a function of p and u. If we indicate the l.h.s. of Eq. fl2.11|) by 
F(x, y, p), the equation for 71 has the form: 



F{x(ji),y[l2(p,7uu)], p} = 0. 



Eq. ( p.5|) then becomes: 



B{p,u) 



C(p,u) 



' d 2 u 
dp\ 



(2.12) 



(2.13) 



where B(p,u) and C(p,u) are given by: 



B{p,u) 
C(p,u) 



V d l2 

;i-0 2 du 



df dF dx df dF dy 



P 



8^2 dx dji 
dF dy d r )2 dF dx 
dy d^2 dji dx dji 



971 dy dj 2 



(2.14) 
(2.15) 



where all the partial derivatives are performed at constant density p. 

We note that in Eq. ( |2.1| ) there is no hard-sphere contribution to the direct correlation 
function outside the repulsive core. As a consequence, the treatment of thermodynamics and 
correlations of the hard-sphere gas that come out of Eq. ( |2.1| ) in the high-temperature limit 
coincide with that of the Percus-Yevick (PY) integral equation. An improved description 
of the hard-sphere thermodynamics at the level of the well known Carnahan-Starling (CS) 
equation of state is desirable in this context, as the slight inaccuracy in the PY treatment 
of the hard-sphere gas affects not only the high-density behavior of the system in study, 
but also the location of its critical point and of the coexistence curve. A non-vanishing 
contribution to c(r) for r > 1 due to hard-core part of the interaction could be taken into 
account by replacing the expression for c(r) of Eq. (|2.1|) with 



c(r) = c H s(r) + K {P, P)w(r) 



r > 1 



(2.16) 



where the hard-sphere direct correlation function chs(?") can be determined, for instance, 
by the Waisman parameterisation |T3|]. This is what has been done in Ref. for the HCY 
fluid. However, in the present case of a two- Yukawa tail potential, this would require one to 
deal with a direct correlation function of three- Yukawa form for r > 1. Although such an 
extension of our treatment appears to be feasible , in order to minimize the complexity of 
the computation, we have not pursued it here. Instead, we assumed for cns( r ) outside the 
core a Yukawa form whose range coincides with that of the repulsive contribution to the tail 
potential, and whose density-dependent amplitude H is set so as to give CS thermodynamics 
in the high-temperature limit. This procedure does not aim at the most accurate description 
of the hard-core contribution to the correlations possible in the context of the SCOZA. But 
it has the advantage of taking into account the hard-sphere thermodynamics beyond the PY 
level without going beyond the two- Yukawa form of c(r) . An analogous procedure involving 
only a single Yukawa term was already implemented for the HCY fluid fllS], where it was 



5 



found to give results very similar to those of the more reliable closure ( |2.16|) both for the 



thermodynamics and the phase diagram. According to this prescription, for r > 1, c(r) has 
the form: 



c(r) 



(H + KAA r . ,. KA 2 r , 
exp[— z\{r — 1)J exp|— Zi\r — 1)\ r > 1 , (2-17) 



while for r < 1 the core condition g{r) =0 holds as before. Here K is the unknown, state- 
dependent amplitude of Eq. fl2.1|) , which vanishes in the high-temperature limit, while if is a 
known function of the density. As the inverse range of the hard-sphere contribution to c(r) is 
locked to zi, the treatment of the hard-sphere gas that comes out of Eq. (|2.17| ) clearly lacks 
the virial-compressibility consistency of the Waisman parameterisation. However, this does 
not affect the consistency between internal energy and compressibility route upon which 
SCOZA hinges. The manipulations that lead to the PDE ( |2.13| ) remain unchanged, the only 
difference being that Eq. Q2.11|) will contain an extra term related to H, which will affect 
the partial derivatives dF/dx, dF/dy that appear in Eqs. ( |2.14| ), ( |2.15| ). 



The PDE ( p,13| ) has been integrated numerically. The initial condition at f3 — and 



the boundary conditions at the ends of the density interval are the same as in Ref. || and 
will not be detailed here. Once u has been obtained by solving Eq. fl2.13| ), integration of u 
with respect to (3 yields the Helmholtz free energy and hence all the other thermodynamic 
quantities. 



III. SIMULATION PROCEDURE 

The principal aspects of the simulation and finite-size scaling techniques employed in 
this work have previously been detailed elsewhere in the context of a similar study of the 
Lennard- Jones fluid. Accordingly we restrict ourself to a brief summary of the methodology 
and refer the reader to reference |16| for a fuller account. 

Grand canonical Monte-Carlo (MC) simulations were performed for the HC2YF model 
of Eq. ( |1 . 1|) . The algorithm used had a Metropolis form |T7j and comprised only particle 



transfer (insertion and deletion) steps, leaving particle moves to be performed implicitly as 
a result of repeated transfers. The potential was cut at a radius r c = 3. Oct, and a standard 
correction term was applied to the internal energy to compensate for the trunction [|17]]. To 
simplify identification of particle interactions the periodic simulation space of volume L 3 
was partitioned into m 3 cubic cells, each of side the cutoff r c . This strategy ensures that 
interactions emanating from particles in a given cell extend at most to particles in the 26 
neighbouring cells. 

System sizes having m = 3,4, 5, 6, 7 and 8 were studied, corresponding (at coexistence) to 
average particle numbers of approximately 230, 540, 1050, 1750, 2900 and 4500 respectively. 
For the m = 3,4,5 and 6 system sizes, equilibration periods of 10 5 Monte Carlo transfer 
attempts per cell (MCS) were utilised, while for the m = 7 and m = 8 system sizes up to 
2 x 10 6 MCS were employed. Sampling frequencies ranged from 20 MCS for the m = 3 
system to 250 MCS for the m = 8 system. The total length of the production runs was also 
dependent upon the system size. For the m = 3 system size, 1 x 10 7 MCS were employed, 
while for the m = 8 system, runs of up to 1 x 10 8 MCS were necessary. 
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In the course of the simulations, the observables recorded were the particle number 
density p = N/V and the energy density u = E/V. The joint distribution p L (p,u) was 
accumulated in the form of a histogram. In accordance with convention ||17|| , we express all 
thermodynamic quantities in reduced units: p* = pa 3 , u* = ua 3 /e, T* = k B T/e. 

Efficient exploration of the phase space was facilitated through use of the histogram 
reweighting technique fllSf . This method allows histogram accumulated at one set of model 
parameters to be reweighted to provide estimates appropriate to another set of not-too- 
distant model parameters. Use of the method permits large areas of phase space to be 
mapped using only a few simulations performed at strategic state points. To facilitate study 
of the subcritical coexistence region, the multicanonical preweighting technique |19j was 
employed. This method employs a biased sampling technique to overcome the free energy 
barrier separating the coexisting phases and thus allow both to be sampled in a single 
simulation run. When combined with histogram reweighting in the manner described in ref. 
rj| , multicanonical preweighting permits an extremely efficient accumulation of coexistence 



curve data. 

The task of estimating the critical point parameters of the model was performed using 
finite-size scaling techniques ||16|| . In brief, the strategy is to match the measured "ordering 
operator distribution" to an independently known universal critical point form appropriate to 
the Ising universality class. The ordering operator itself is defined as M. oc (p* + su*), where 
s is a non-universal "field mixing" parameter, which is finite in the absence of particle-hole 



symmetry, and which is chosen to ensure that p(Ai) is symmetric in Ai ||20|| . For sufficiently 



large L, the matching should occur at the critical point parameters. For small L, however, the 
apparent critical temperature obtained by this matching procedure is subject to systematic 
errors associated with corrections to finite-size scaling. To deal with this, we extrapolate 
to the thermodynamic limit using the known scaling properties of the corrections, which 
are expected to diminish (for sufficiently large system sizes) like L~ 9 / u |L6 |, where 8 is the 
correction to scaling exponent and v is the correlation length exponent. The extrapolation 
has been performed using a least squares fit to the data for the four largest system sizes. 
The results of the extrapolation are shown in Fig. |2|, from which we estimate T* = 1.295(10). 
The associated estimate for the critical density is p* = 0.310(1) and for the reduced chemical 
potential is p* = -3.588(30). 

In addition to the phase coexistence data, we have also measured the form of the radial 
distribution function g(r) for a number of state points, corresponding to reduced tempera- 
tures T* = 2, 1.5, 1, and reduced densities p* = 0.4, 0.6, 0.8. In the following Section, the 
simulation results for the coexistence curve, critical point properties and radial distribution 
function of the HC2YF are compared with the theoretical predictions. 



IV. RESULTS AND DISCUSSION 

The SCOZA coexistence curve has been determined by equating the pressure P and the 
chemical potential p on the low- and high-density branch of the subcritical isotherms. As 
we specified above, these quantities were obtained via the energy route by integrating u 
with respect to (3. The advantage of doing so is that one does not have to circumvent the 
forbidden region bounded by the spinodal curve, i.e. the locus of diverging compressibility, 
in order to obtain P and p on the high-density branch. On the other hand, because of the 
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compressibility-internal energy consistency of the theory, this is fully equivalent to using 
a mixed path combining integration of the inverse compressibility 1/Xred with respect to p 
and integration of u with respect to (3, as is often done in calculations based on integral 
equations [@. The coexistence curve in the density-temperature and in the temperature- 
chemical potential plane is shown in Figs. |3] and |] respectively. Reduced units have been 
used throughout. The SCOZA results obtained both by the closure (|2.1| ) and by the mod- 
ified version ( |2.17| ) are compared with the MC data obtained in this work. We have also 
plotted the coexistence curve predicted by HRT and by the energy route of the lowest- 
order gamma-ordered approximation (LOGA) |HJ, also known as optimized random-phase 
approximation (ORPA) |22[, together with the spinodal curve given by LOGA/ORPA com- 
pressibility route. It appears that the SCOZA yields the most accurate determination of the 
coexistence curve among the theories considered here, provided Eq. Q2.17| ) is used in order 
to describe hard-sphere thermodynamics at the CS level. The predictions for the critical 
point are compared in Tab. |. The error on the SCOZA critical density and temperature is 
below 1%. Actually, the critical temperature predicted by SCOZA with closure ( |2.1| ) and 
PY hard-sphere thermodynamics is even closer to the MC result, but this is most likely to 
be accidental since, as stated in Sec. II, on the basis of our previous calculations on the 
HCY potential [Qj we expect that a treatment of the HC2YF based on closure ( |2.1(j| ) and a 
c(r) of three- Yukawa form will in fact deliver results almost identical to those found here by 
Eq. (ETT7D . We note that, like the HRT or the energy-route LOGA/ORPA, the SCOZA gives 
a coexistence curve that goes right up to the critical point. This feature is not shared by 
other integral-equation approaches, such as the MHNC or the HMSA theories [0]. Both the 
SCOZA below the critical temperature and the HRT yield non-classical critical exponents; 
for the exponent /3 coex which describes the curvature of the coexistence curve the SCOZA 
gives /3 coex = 7/20 = 0.35 G3J, while according the HRT /3 coex ~ 0.345 [O, the best theoret- 



ical estimate being (3 coc ^ ~ 0.327 ||24|| . Moreover, in the SCOZA, as well as in the HRT, the 
critical point, identified as the top of the coexistence curve, coincides with that determined 
by locating the divergence of the isothermal compressibility as given by the compressibility 
sum rule. As a consequence, at the critical point one has both coalescence of the vapor and 
liquid phases and occurrence of long-range correlations, as expected. Because of the lack of 
thermodynamic consistency, this is not the case with LOGA/ORPA, in which the critical 
point obtained by the compressibility route corresponds to the top of the spinodal curve 
shown in Fig. |3|. 

We now consider the radial distribution function g(r). The SCOZA g(r) has been deter- 
mined using both closure (Q and (pTQ . In Fig. | the SCOZA and LOGA/ORPA results 
for two different states are compared with MC simulations. We see that, while the overall 
agreement is satisfactory, the behavior near contact is not reproduced very well. It is unlikely 
that this can be traced back to the slight inaccuracy in the treatment of the hard-sphere 
contribution to c(r) for r > 1 entailed by Eq. ( |2 . 1[ ) or ( |2.17| ). In fact, such a contribution 
is accurately taken into account in LOGA/ORPA, which nevertheless does not appear to 
perform better than SCOZA. In particular, using closure ( |2.17| ) in SCOZA produces results 
for g(r) which are undistinguishable from those of LOGA/ORPA, at least for the states we 
investigated. Moreover, we checked that for the densities considered here the hard-sphere 
g(r) obtained from closure ( [2. 17[ ) in the high-temperature limit is nearly superimposed to 
that given by the Waisman parameterisation. On the other hand, both in SCOZA and in 
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LOGA/ORPA the contribution to c(r) due to the tail potential w(r) is bound to follow 
the profile of w(r) itself for all r's. Such a form is best suited for interactions that are 
slowly varying on a lengthscale of the order of the particle size, while the two- Yukawa tail 
of Eq. (|1 . 1| ) changes quite steeply near its minimum. As is well known, such a problem is 
usually dealt with by splitting the potential according to the WCA prescription [[□]], but 
implementing this procedure in the SCOZA seems very artificial, as it would hinder the ana- 
lytical tractability of the system, which is the very reason why the HC2YF interaction ( |1 . 1| ) 
was considered in the first place instead of the LJ potential. 

Since on the basis of the comparison with MC results the SCOZA proves to accurately 
reproduce the coexistence curve of the HC2YF, it could be worthwhile using it to assess 
the ability of different HC2YF parameterisations of the LJ potential to reproduce the LJ 
coexistence curve. To this end, in Fig. ^| we have compared the coexistence curve of the LJ 
potential given by MC simulations [25] with the SCOZA results using closure ( 2.17|) for the 
parameterisations proposed by Sun 0,[7|] and by Foiles and Ashcroft |J, together with that 
of Eq. ALT] ). We see that the HC2YF form of Eq. (Q, closely resembling that of Ref. [|, 
gives the best representation of the LJ potential among those considered here, at least as 
far as the liquid- vapor phase diagram is concerned. 

In summary, we have presented SCOZA results for the phase diagram and the correlations 
of a two- Yukawa parameterisation of the LJ potential. These have been compared with the 
predictions of other theories and with new MC simulation data. The comparison shows that 
the SCOZA provides a very accurate coexistence curve, with a critical point that differs 
from the prediction of MC supplemented by finite-size scaling techniques by less than 1%. 
At the same time, the two- Yukawa potential we studied, closely resembling that previously 
proposed by Kalyuzhnyi and Cummings 0, reproduces satisfactorily the LJ coexistence 
curve. In view of the very good performance of SCOZA in predicting the liquid-vapor phase 
diagram of simple fluids, we think that it could be worthwhile to generalize the solution 
procedure, so that it will not be confined anymore to tail potentials of Yukawa form. This 
could also allow one to implement SCOZA with a more sophisticated closure than the forms 
and (|2~T7|) used here. 
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TABLES 





MC 1 " 


SCOZA* 


SCOZA* 


HRT 


LOGA| n 


LOGA^ omp 


Pi 


0.310(1) 


0.307 


0.304 


0.310 


0.314 


0.328 


± C 


1.295(1) 


1.304 


1.293 


1.316 


1.352 


1.071 



TABLE I. Critical density and temperature (in reduced units) for the HC2YF. f : MC simula- 
tion performed in this work. SCOZA result using 2- Yukawa c(r) of Eq. (|2,17|). •: SCOZA result 



using 2- Yukawa c(r) of Eq. (|2T|) . \: LOGA/ORPA result using energy route, o: LOGA/ORPA 



result using compressibility route. 
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FIGURES 




FIG. 1. Pair potentials of the interparticle interaction given by the LJ fluid (dotted curve) and 



the HC2YF (full curve) of Eq. (1.1) with amplitude and range parameters specified in the text. 
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1.294 



H 1.292 



1.29 



1.288 








0.01 



0.02 

-(l+6)/v 



0.03 



0.04 



FIG. 2. The apparent reduced critical temperature (as defined by the matching condition de- 
scribed in the text and in [jlfj), plotted function of IrC^ 1 )/", with 9 = 0.54 and v = 0.629. 
The extrapolation of the least squares fit to infinite volume yields the estimate T* = 1.295(1). 
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FIG. 3. Coexistence curve of the hard-core two- Yukawa fluid in the density-temperature plane. 



Density and temperature are in reduced units. Solid line: SCOZA with closure ( 2.17 ) and Carna- 
han-Starling reference system thermodynamics (see text). Dashed line: SCOZA with closure ( |2.1| ) 
and Percus-Yevick reference system thermodynamics. Dotted line: LOGA/ORPA (energy route). 
Circles: HRT. Squares: MC simulation performed in this work. Cross: MC critical point. 
Dot-dashed line: LOGA/ORPA spinodal curve (compressibility route). 




FIG. 4. Coexistence curve in the T*-fj,* plane. Symbols as in Fig. 
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FIG. 5. Radial distribution function of the HC2YF at T* = 1.5, p* = 0.4 (upper panel) and 
T* = 1, p* = 0.8 (lower panel). Solid line: SCOZA with closure fl2iP . Crosses: LOGA/ORPA. 
Circles: MC simulation. The result of SCOZA with closure ( |2.17 ) are undistinguishable from those 
of LOGA/ORPA. 
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FIG. 6. Comparison between the LJ and the HC2YF coexistence curve for different HC2YF 
parameterisations. Lines: SCOZA results for the HC2YF of Eq. fll.lD with amplitude and range 
parameters by Kalyuzhnyi and Cummings (solid line), by Foiles and Ashcroft M (dotted line), 
and by Sun BLR] (dashed line). Triangles: MC simulation results for the LJ coexistence curve pF 



Diamond: MC simulation result for the LJ critical point [26|. 
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